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^ Abstract 

T— I In the present paper, based on the notion of niarginaUy stable modes of vector 

^_^ perturbations in the AdS soliton back ground we investigate the magnetic field effect 

^S^ on holographic p- wave insulator /superconductor transition in the probe limit. We 

perform explicit analytic calculations considering both Schwarzschild AdS soliton 
and Gauss Bonnet AdS soliton back grounds and obtain a unique relation between 

4^ the chemical potential and the magnetic field near the critical point. We also extend 

our analysis for p + ip- wave back ground in the presence of external magnetic field. 

cn In both the cases it is observed that the non abelian model exhibit droplet solutions. 

^ Moreover it is also found that the increase in the value of the external magnetic 

r-^ field essentially reduces the size of the droplet and thereby makes the condensation 

^— I harder. 

^ 1 Introduction 

Kj) For the past few years AdS/CFT duality [l]-[3] has received significant attention due to its 
'O several remarkable applications in various areas of physics particularly towards the under- 
^ standing of several strongly coupled phenomena in the usual condensed matter systems. 
^ Among its various remarkable achievements, the gauge/gravity duality has successfully 
explained various features of so called high Tc superconductivity in the usual condensed 
matter physics. These are known as holographic superconductors |l]-[6]. The main idea 
behind such theories lies in the fact that if the temperature of the black hole in the AdS 
space falls below certain critical temperature (T < T^) the black hole would become un- 
stable to develop a nontrivial hair profile through the mechanism of abelian symmetry 
breaking in the bulk gravity theory. This essentially triggers superconductivity/super- 
fluidity in the boundary CFT through the mechanism of spontaneous breaking of a global 
(or sometimes weakly gauged) U{1) symmetry. For T > T^ the black hole possess no hair 
and the corresponding gravity theory stands for the dual of a conductor. Therefore these 
models essentially describe conductor/superconductor transition in the boundary CFT 
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Besides having the usual description of conductor/superconductor phase transition, 
one may also generate so called insulator/superconductor phase transition in the boundary 
CFT through the appropriate modifications of the bulk gravity theory, i.e; considering 
the gravitational theory in an AdS soliton back ground |15]. The AdS soliton that has 
been first investigated in [36] could be obtained through a double wick rotation of the 
usual AdS Schwarzschild solution. In the language of gauge/gravity duality the AdS 
soliton is dual to a confined field theory with a mass gap that essentially corresponds to 
an insulating phase ^\. In |15], for the first time, it had been argued that in the presence 
of a chemical potential (/i) in the AdS soliton back ground, apart from having the usual 
first order insulator/conductor phase transition between AdS soliton and AdS black hole 
[18] . one may have an additional (second order) phase structure through the appropriate 
tuning of /i. It has been observed that for /i > fic, where fj,c is some critical value of the 
chemical potential, the AdS soliton back ground becomes unstable to develop a non trivial 
profile of hair which could be interpreted as insulator/superconductor phase transition in 
the language of boundary field theory. This is kind of quantum phase transition where the 
chemical potential (/x) plays the role of doping (x) in the ordinary material, like cu-prates. 

Since the discovery of such an exciting phase structure, till date a number of inves- 
tigations have been performed in various directions. In |39],[5nj the analysis was further 
extended taking into account the back reaction of the matter fields. It was observed that 
one may generate a completely new type of phase structure depending on the strength of 
the back reaction. In order to back up the previous numerical results [51], full analytic 
calculations were carried out in [5^ , [H3] both for the Scwarzschild AdS and Gauss-Bonnet 
AdS back grounds. In [51] the authors studied the insulator/superconductor phase tran- 
sition via Stuckelberg mechanism and unveiled reach physics in various phase transitions. 
Response to a Wilson line on the circle and to a magnetic field perpendicular to the 
non compact directions have been explored in [55] where the authors discover Aharonov- 
Bohm like effects during insulator/superconductor transition. A holographic model of SIS 
Josephson junction have been successfully constructed in [5S] and the dependence of the 
maximal current on the width of the junction have been investigated. Based on the no- 
tion of marginally stable modes of scalar /vector perturbations in the AdS Schwarzschild 
background, the authors in [57], [58] find that these modes can actually indicate the onset 
of insulator/superconductor transition which is also compatible with earlier observations 
made by Gubser [n],[22]- Apart from these studies, very recently some research papers 
were written on insulator/superconductor transition in the context of non linear electro- 
dynamics as well as entanglement entropy [59]- 162] . 

From the above list of references, it is quite interesting to note that most of the papers 
that are written so far are mostly based on holographic s- wave model. On the other 
hand, a surprisingly lesser effort is given towards the study of non abelian model of in- 
sulator/superconductor transition [63] , particularly the magnetic response in non abelian 
model of insulator /superconductor transition is still unexplored. The present paper aims 
to fill up this gap along with addressing some other major issues from a different perspec- 
tive. Based on analytic calculations, the present work reveals the importance of marginally 
stable modes of vector perturbations (both for the Schwarzschild AdS and Gauss-Bonnet 
AdS soliton back grounds) in order to construct non abelian superconducting droplets 
(solutions that are confined with in a finite region and decays rapidly at large distances) 
during insulator/superconductor transition in presence of an external magnetic field. It 
is noted that the size of the droplet gradually reduces to zero as the strength of the 



magnetic field is increased. Furthermore, it is found tfiat near tfie critical point of phase 
transition there exists a simple algebraic relation between the chemical potential (/i) and 
the magnetic field [B) which indeed reveals the fact that the increase in the strength of 
the magnetic field makes the condensation harder to form. On top of it, for Gauss-Bonnet 
gravity we find this simple relation to be dependent on the Gauss-Bonnet coupling (a) 
which is also an interesting observation in itself. 

It interesting to note that the droplet solutions that are obtained throughout this 
paper is in line with the earlier observations during conductor /superconductor transitions 
where both the vortex and droplet solutions have been constructed for the s- wave model 
[35], [37] -[38]. From the technical point of view the absence of vortex states in the non 
abelian model could be explained as follows: In the non abelian holographic model the 
Maxwell field enters as a U{1) subgroup of the full SU{2) symmetry group. Therefore 
it is not just simply coupled with the condensed charged fields through gauge covariant 
couplings like in the s or d- wave case. Thus technically speaking the nature of this 
difference lies in the nature of the minimal coupling. For real life materials the high Tc 
type II superconductivity is believed to be generated through a d- wave pairing mechanism 
which results into vortex states when magnetic field is applied externally. There are also 
real world p + ip-wave models that exhibit vortex solutions but that are clearly different 
from the holographic p + ip- wave model studied in this paper. On the other hand, 
one could interpret our model as a holographic model of super-fluidity and try to find 
connections with real life super-fluid droplets [M]. 

Before going further, let us briefly mention about the organization of the paper. In 
section 2, we discuss the importance of marginally stable modes of vector perturbations 
in order to study non abelian superconducting droplets. The calculations are carried out 
both for the Schwarzschild AdS and Gauss Bonnet AdS solitonic back grounds. In section 
3, the droplet solutions have been constructed for the p + ip- wave back ground. Finally, 
the paper is concluded in section 4. 

2 Droplet solution via marginally stable modes 

By marginally stable modes in a (black hole) space time we mean perturbations in space 
time that do not back react on other fields and depend only radial coordinates (r). In other 
words, marginally stable modes are essentially the quasi normal modes whose frequen- 
cies go to zero {u = 0) near the critical point of the phase transition. The importance of 
marginally stable perturbations in case of holographic model of conductor /superconductor 
transition has been first discussed by S. S. Gubser in his pioneering papers [9j,[22j, where 
he discussed both the abelian and the non abelian model of superconductors. Later 
on in [5S] the authors investigate the role of marginally stable modes in AdS soliton 
back grounds and find that these modes indeed indicate the onset of so called insula- 
tor/superconductor transition near the critical point. 

Motivated from the above studies in the subsequent analysis we first investigate the 
role of marginally stable modes of vector perturbations for the p- wave back ground 
(where we consider perturbations only in one spatial direction r^dx) in presence of an 
external magnetic field and in the next section we generalize our calculations considering 
perturbations both in the r^dx and r'^dy directions. 



2.1 Schwarzschild AdS soliton 

In this section we aim to construct holographic droplets for the non abelian case based on 
the idea of marginally stable modes. In order to do that we first note that the five dimen- 
sional SU{2) Einstein- Yang-Mills action in presence of a negative cosmological constant 
may be written as |63j, 



S = / d^x^/^ 



-iR-X)--F''F^^"' 



(1) 



where F^,y is the field strength of the non abelian gauge field. 

Before proceed further let us first note that through out this paper calculations are 
performed in the probe limit and we treat marginally stable perturbations as a probe into 
the neutral AdS soliton which may be expressed a.¥\ j46] , 

ds^ = -$-- + r\~de + dp" + p'dO^) + r^g{r)dx^ (2) 

where, g{r) = 1 — ^ and L = 1. 

This solution can be obtained from a five dimensional AdS black hole solution followed 
by two successive wick rotations. The boundary CFT that we are studying actually lies 
in a space which has the topology i?^'^ x S^. 

In order to proceed further we choose the following ansatz for the gauge sector close 
to the critical point of the phase transition line (i,e; /z ~ /Xc and ^ ~ 0), 

A = T^ifijt + -Bp^de) + ^{t, r, X, pydp (3) 

where r^ could be identified as the generator of the electromagnetic U{1) subgroup of 
the full SU{2) gauge group. Here _B is a constant magnetic field which is related to the 
electromagnetic vector potential. 

With the above choice, the equation of motion for if) turns out to be. 



dr{r'^g{r)dri)) - -d^ip + ::77:x^x^ + ^'^- -^^ = 0- W 



1 1 //2 R2n2 

In order to solve the above equation Q, we choose the ansatz of the following form 



^ = F{r,t)H{x)Uip) = F{r,t)H{x)^ (5) 



U{p) = p=-=. (6) 

V -D 



where we have actually set. 



From the above analysis it is indeed clear that the solutions we are looking for is 
actually confined with in a finite circular region whose radius p oc -^. Thus for sufficiently 
large magnetic field the model essentially describes a superconducting droplet for which 
the condensate is non zero only at the origin where the core of the droplet is located 

mm-m- 



^We are working in the polar coordinates (p, 9). 



Substituting the above ansatz ([s]) into ^ we find the following set of equations, 

"■^^""^ dr{r'g{r)drF{r,t)) - .^P-d^F{r,t) + g{rM - B) = k" (7) 



F(r,t) ^ "^ ' ^ ' " F{r,t) 
and, 

dlH = -K^H. (8) 

The solution of ([s]) may be noted as, 

H{x) = expiiKx) (9) 

where, we identifjjjfi; = 2n from the periodicity propertjjjof H{x)- 
Finally, equation ([T]) turns out to be, 

d'rF{r,t) +(- + ^-P-) drFir,t) - -^^d^F{r,t) + ^ ( f^l - B - ^) F{r,t) = 0. 

(10) 

As a next step, we substitute F{r,t) = exp{—iwt)R(r), which essentially reduces ( |lO| 
into a radial equation in the variable R{r), 

,2nr^^ , /'3 , 9\r)\^ ^,, , / , , , ^ 4^^ \ i?(r) 



dR{r) +[- + '-ri 9rR{r) + u;^ + /i - 5 - -- ^y- = 0. (11) 

Since in the present analysis we are interested in perturbations that are marginally 
stable, therefore we set a; = 0. Also we redefine our variable as 2; = 1/r. With this 



redefinition equation (11) turns out to be. 



,2o/-> , /'s'W l\a u,.^ , I. .2 o 4n2\ fl(.-) 



As a next step, we introduce a trial function VL{z) as, 

i?(z)|,_0~(O2)^'fi(^) (13) 

with the boundary conditions for VL{z) as ^2(0) = 1 and Vt (0) = 0. 



Substituting (13) into (12) the equation of motion for Vt[z) turns out to be 

dlVL{z) + p{z)d,n{z) + q{z)n{z) + \^w{z)Vl{z) = (14) 

where we have set X^ = fj,^ — B along with, 

3 — 7z'^ 
p{z) 

q{z) 



z{\-z^) 
-8z^ 4^2 



1 - z^ g\z) 
w{z) = ^^. (15) 



^We have actually an infinite tower of marginally stable modes corresponding to diff'erent values of n. 
It is generally expected that the lowest mode of excitation {n — 0) will be the first to condense and all 
other modes do not matter to the phase transition. 

3iJ(x) - i/(x + tt) 



It is now easy to convert (14) into the standard Strum-Liouville (SL) eigen value 
equation [14j, 

d,{T{z)n'{z)) - Qiz)n{z) + x^p{z)n{z) = o (i6) 

with the following identifications, 

T{z) - 

Q{z) - 

P{z) = z\ (17) 

The corresponding eigenvalue (A^) is obtained through minimizing the following func- 
tional, 



.3(1- 

z^{8z 




z^)+An') 


^3 


1- 


-z* 



\'[n{z)] = ^,i-B 



jldz{T{z){^{z)f + Q{z)n\z)) 
CdzP{z)n\z) 



In order to estimate the R.H.S. of (18) we take VL{z) = 1 — az"^, which finally yields, 

III- B = 5.13933 (19) 

with a = 0.338114. 

This is an important result which states that the chemical potential (/x ~ fic) and the 
magnetic field {B ~ Be) may be connected to each other via a simple algebraic relation 
near the critical point of the p- wave insulator/superconductor phase transition. Because 
of the presence of the negative sign it is easy to note that the increase in the value of the 
external magnetic field essentially increases Hc which results in a harder condensation. 

In the next section we wish to verify this relation for the Gauss-Bonnet AdS soliton 
back ground and we will try to examine the effect of so called higher derivative coupling 
on the non abelian condensate in presence of external magnetic field. 

2.2 Gauss-Bonnet AdS soliton 

It has been noted earlier in [53], that the Gauss-Bonnet coupling (a) indeed affects the 
insulator/superconductor transition. Increase in the value of a essentially increases fie 
which results in a harder condensation. In the present section we aim to investigate 
the situation in presence of an external magnetic field. The motivation of the present 
calculation is to see how the simple algebraic relation that is obtained in the previous 
section now depends on the value of the coupling a. 

We start with the five dimensional AdS soliton in the Gauss-Bonnet gravity which (in 
polar coordinates) takes the following form [53j|, 



ds' = -4^^ + \hdt' + dp' + p'de') + ^-^dx' (20) 

z'^g[z) z^ z^ 



where, g[z) = ^ 1 - J\ - 4a(l 



with a as Gauss-Bonnet coupling. 



The asymptotic form of the metric is given by. 



^WI.^o~^[l-v/r^4^ (21) 



from which we identify the effective AdS length to be, 

2a 



'"" 1 - yr^4^' 

We take the same ansatz (Is]) for the gauge field (A^) and set, 

^(^, z, X, p) = Riz)S{x, p)exp{-iut). 



(22) 



(23) 



We proceed in the same way as we did in the previous section. Since we are interested 
in marginally stable modes therefore we will be interested in the equation of motion for 
R{z). Finally we set w = 0, which yields, 



where, $,{z) 



dlR{z) + ^{z)d,R{z) + 



l-Aa-Aaz^~^l-Aa(l-z^) 
z |"^l-4a(l-24)_i+4Q,(i_^4)l 



pI-b~ 



An 



2 n 



9[z) 



9{.z) 







(24) 



Our next step would be to convert (24) into the standard SL eigen value equation. In 
order to do that we first set 

R{z)l^,^{02)z^n{z) 

which yields. 



dln{z) + p{z)d,n{z) + q{z)n{z) + yw{z)n{z) = 



(25) 
(26) 



where we note, 



p{z) 



q{z) 



3. 


/I - 4a(l - z^)-3 + 12a - 20az^ 






z 


^1 - 4a(l -z^)-l + 4a(l - z^) 



■IQaz'^ 



l-4a(l-z4) -l + 4a(l - ^^ 



4n2 
g^{z) 



^iz) = g (z) 



(27) 



with, X"^ = nl- B. 



This equation (26) may be converted into a standard SL eigen value equation, 

d,{T{z)Q'{z)) - Q{z)Q{z) + X^P{z)Q{z) = 



(2^ 



with the following identifications, 

r3 



Z" 

2a 



1 - 4afl - z^ 



Tiz) 

z.a L 

Q{z) = -T(z)g(z) = 8z^(l + 2a(l-/)) + C(a2) 
P{z) = T{z)w{z) = z\ 



z%l - z^){l + a{l - z^)) + 0{a' 



(29) 



Since the value of the higher derivative coupling (a) is very small therefore for the present 
analysis it will be sufficient to consider the leading possible effect of the higher derivative 
correction to the onset of p- wave insulator/superconductor transition. 
The eigen value could be obtained by varying the following functional. 



A [Q{z)] =fi^-B = -^ = ({a, n) (30) 

Jo dzP{z)Q^{z) 

In order to estimate the above functional integral we consider the trial function to be 
of the form Q{z) = 1 — az"^, where a is determined through the minima of the functional. 
In the following (Table 1) we tabulate various values of ({a, n = 0) corresponding to 
different choice of the coupling parameter (a). 

Table 1: Values of ( for different choices of a 



a 


a 


CSL (= /i^ - B) 


0.0001 


0.33809 


5.13978 


0.1 


0.314773 


5.59003 


0.2 


0.29277 


6.03793 



From the above table (1) it is indeed clear that in the presence of higher derivative 
corrections to the bulk action, the simple algebraic connection between the chemical 
potential (/Xc) and the magnetic field (B) is now dependent on a. Furthermore we note 
that given a fixed value of B, the corresponding value of fie increases as we increase 
the value of the Gauss-Bonnet coupling (a). This indeed suggests the onset of a harder 
condensation and which is expected from the earlier observations 



3 Droplet solutions in p + ip- wave back ground 

In the present section, we would like to investigate the magnetic response in holographic 
p + ip- wave superconductors considering the probe limit. In the new coordinate z = 1/r, 
the metric of the space time may be written as, 

ds^ = -p-- + ^{-dt" + dx^ + dy^) + ^rfx' (31) 

z'^g[z) z^ z^ 

where, the explicit form of fy'(z) depends on the particular choice of space time. 

In this fixed back ground, the corresponding equation of motion for the gauge field 
turns out to be, 

--d^{^^F''^'')^e"'^^A\F'^^'' = ^ (32) 



with the choice e"^^ = 1. 

In order to proceed further we choose the following ansatz for the gauge field in the 
p + ip- wave back ground as, 

A = (0(x, y, z)dt + A^(x, y, z)dy)T^ + w{x^ y, z)(T'^dx + r'^dy) (33) 

where the condensate of w bosons essentially breaks the original rotational invariance of 
the system and thereby triggering superconductivity in the boundary theory. 

Substituting the ansatz (33) into (32) the following set of equations for 0, A^ and w 



may be respectively obtainedas 



a , I 9[z) 1\ -. , , 1 /a2 , a2N. 2w' 



d.r+{^-fi-- 5.0+ —-(9,^ + 9,^)0- —f = (34) 

\g{z) z) g{z) y g{z) 

8 



and. 



d. ( ^5.4 ) + -dlAl + —d^w - -w'Al = 

/C / /C iC A^ 






(35) 



(36) 



From the structure of these equations it seems to be quite difficult to have exact ana- 
lytic solutions for these equations. However, one can always make perturbative expansion 
of the fields near the critical value of the magnetic field strength [B ~ Be) where there is 
no condensation (w = 0), 



A, 



y 
w 



0(0) +^0(1) + .... 

Af)+eAf) + .... 



Bc-B 
Be ' 



where the perturbation expansion is carried out on the parameter e 

slightly below B^.. Here (j)^'^^ and Ay are zeroth order solutions which we set, 

= /ic and, A^y^^^ = B^x 



(37) 



when B is 



(3^ 



in order to make the solutions consistent with the Newmann type boundary conditions 
for the fields [52j. 

Substituting (37) into (36), and neglecting term^~ 0{£w) we obtain. 



2zd, (^d,wi\ + {dl - By - 3B,)w, + d'yWi + 2/i>i = 



In order to solve (39), we choose the following ansatz, 

wi{x,y,z) = X{x)Y{y)Z{z). 



Finally, substituting (40) into ( p9| ), we arrive at the following set of equations, 

{-dl + x^ + 3-k^)X{x) 







and, 



d!Z{z) + 



g{z) zj V ^ V ^ / g(^z) 







(39) 

(40) 

(41) 
(42) 

(43) 



wherdj, we set x = \/BcX and y 



^cV. 



Solutions for (41) and (42) may be expressed as, 

Y ~ exp{ipy) 



(44) 



Near the critical point, w ^ 0, therefore we may ignore the quadratic terms. 
^We choose the constants fc^ and p^ such that k'^ + p^ = 2. 



and, 

X{x) = e-^Hnix) (45) 

respectively, where Hn{x) are standard Hermite functions. 

Before we proceed further, let us pause for some time and try to extract physics from 



the above results. First of all, and most importantly, from (45) we note that the solution 
in the x direction does not depend on p. This situation is completely different from the 
s- wave case, where the solutions obtained for x directions essentially depend on p, which 
results in a vortex lattice solutions with a periodicity in the x direction [35]. Thus we 
see that for the non abelian insulator/superconductor transition we can not have vortex 
lattice solutions like we have in the abelian case [S3]. The reason for this is that in the 
non abelian case we have a different scenario where the the Maxwell sector appears as a 
U{1) subgroup of the full SU{2) symmetry group and essentially is not minimally coupled 
with the condensed charged fields like in the abelian case. On top of it, we see that for 



large values of the magnetic field, or at large distances the solution (45) rapidly dies out. 
This on the other hand suggests that these are droplet solutions. Finally, from (43) we 
note that, we may easily convert it into a standard SL form that we did earlier and as a 
consequence of that one can easily show fil — Be = constant. 

4 Summary and final remarks 

In the present paper, based on the notion of marginally stable modes of vector perturba- 
tions in an AdS soliton back ground, several properties of non abelian model holographic 
insulator/superconductor transition have been investigated in presence of an external 
magnetic field in the probe limit. A number of interesting results have been obtained in 
this regard. First of all we have found that the non abelian model exhibit only droplet 
solutions. One of the (technical) reason for this is that in the non abelian case we do not 
encounter a minimal coupling between the Maxwell field and the charged condensed field 
like in the abelian model. As a result of this, the superposition of droplet solutions do not 
give rise to any vortex lattice solution. Similar situations were also encountered earlier in 



On top of it, we have found a simple mathematical relation between the chemical 
potential and the external magnetic field near the critical temperature. This relation in 
itself is very interesting since it implies that an increase in the magnetic field will always 
be accompanied by a corresponding increase in the (critical) chemical potential, which 
in fact suggests the onset of a harder condensation. Moreover, for the Gauss-Bonnet 
gravity we find this relation to be dependent on the Gauss-Bonnet coupling itself. In the 
future one can make further probe into the non abelian model of insulator/superconductor 
transition by studying this model in presence of back reaction and particularly one should 
try to see the effect of back reaction on the holographic condensate in presence of external 
magnetic field. 
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